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(54) Procede d'annulation d'echo et annuleur d'echo mettant en oevre un tel procede 



(57) On applique un signal de reception (x(n)) en en- 
tree d'un systeme generateur d'echo (1 ), on recueille un 
signal d'observation (z(n)) en sortie dudit systeme, on 
filtre le signal de reception pour produire un signal (y(n)) 
de prediction d'echo, et on soustrait le signal de predic- 
tion d'echo du signal d'qbservation pour produire un si- 
gnal d'erreur (e(n)). On decompose le signal d'observa- 
tion (z(n)) en M composantes spectrales se rapportant 
a des fr6quences differentes du spectre, on obtient M 



composantes spectrales correspondantes du signal de 
prediction d'echo (y(n)), on calcule M grandeurs ayant 
chacune une valeur k variation progressive representa- 
tive d'un degre de similitude entre, respectivement, Tune 
des M composantes spectrales du signal d'observation 
et la composante spectrale correspondante du signal de 
prediction d'echo, et on opere au moins une mise en 
forme spectrale du signal d'erreur (e(n)) en fonction des- 
dites M grandeurs calculees pour adapter le filtrage et/ 
ou pour delivrer le signal de sortie de I'annuleur d'echo. 
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Description 

La pr6sente invention conceme le domains de I'annulation d'6cho. Elle concerns notamment, mais non exclusi- 
vement, des annuleurs d'echo acoustique pour installations teiephoniques mains libres. 

5 Un annuleur d'echo, dont la figure 1 montre un schema de principe, traite d'une part un signal de reception x(n) 

provenant d'un disposttif distant et applique k une entree d'un systfeme local 1, et d'autre part un signal d'observation 
z(n) issu du systeme 1 , afin de produire un autre signal e(n) qui est retoum6 vers le dispositif distant par I'intermediaire 
d'un canal de communication. On suppose que les signaux x(n) et z(n) sont des signaux numeriques, echantillonnes 
par exemple k 8 kHz pour de la parole en bande telephonique, ou k 16 kHz en bande eiargie. 

10 L'annuleur d'echo represents sur la figure 1 est destine k une installation tei6phonique mains libres. Le systems 

local 1 comprend ainsi un convertisseur numerique-analogique 2 auquel est applique le signal de reception x(n), un 
haut-parleur 3 recevant la sortie analogique amplifies du convertisseur 2, un microphone 4 qui capts la parols locals 
P mais 6galement du bruit ambiant B et des composantss d'echo EC du signal ds reception duss aux couplages 
acoustiques sntre Is haut-parlsur 3 st le microphons 4, st un convsrtisssur analogique-numeriqus 6 qui produit Is 

is signal d'observation z(n) en num6risant la sortie amplifies du microphone 4. 

L'annuleur d'6cho comprend un filtre num6riqus 7 qui d6livrs un signal ds prediction d'6cho y(n) k partir du signal 

de r6ception x(n). Cs filtrs, qui sst k reponse impulsionnelle finie (RIF) dans le cas considers en reference a la figure 
1 , modeiise la reponse des trajets d'6cho au sein du systems 1 . Le signal cferreur e(n) est d61ivre par un soustracteur 
8 qui retranche le signal de prediction d'6cho y(n) du signal d'observation z(n). 

20 pour tenir compte de la variability des trajets d'echo, l'annuleur d'echo est souvent adaptatif , c'est-^-dire que les 

coefficients du filtre 7 sont regies dynamiquement par un module 9 de controle d'adaptation mettant en oeuvre un 
algorithme d'adaptation approprie. En g6n6ral, I'adaptation des coefficients est effectuee en fonction du signal d'entree 
x(n) et du signal d'erreur e(n). On connatt de nombreux algorithmes d'adaptation possibles, notamment diverses ver- 
sions de I'algorithme du gradient. 

25 Frequemment, l'annuleur d'echo adaptatif est associe k un detecteur d'activite vocals et k un detecteur de double 

parole non reprsssntss sur la figure 1. Ces detecteurs sont utilises pour inhibsr I'adaptation des coefficients du filtrs 
7 lorsqus Is signal de reception x(n) ne comporte aucune activity vocale, et lorsqus dss activites vocalss dscorrslsss 
sont pr6sentes dans les signaux x(n) et z(n). En effet, dans ces deux cas de figure, le signal z(n) ne reprSsente pas 
une bonne observation des trajsts d*6cho exploitable par le module de controls d'adaptation 9. On sait que le detecteur 

30 de double parole est la source de nombreuses difficultes pour le concepteur d'un annuleur d'echo. En particulier, ces 
detecteurs sont frequemment perturb6s lorsqus Is bruit ambiant B sst rslativsmsnt important. En outrs, la discrimination 
des situations de double parole est souvent delicate, et le contr6le en tout-ou-rien de I'adaptation base sur cette dis- 
crimination nuit k la bonne convergence de I'algorithme. 

II existe difterentes formes possibles pour realiser le filtrage du signal de reception x(n), qui sont k leur tour com- 

35 patibles avec difterents types d'algorithmes d'adaptation. Dans le cas classique de la figure 1, on utilise un filtre RIF 
7 realisant la convolution entre le signal de reception et la reponse impulsionnelle estimee, de longueur N echantillons, 
des trajets d'echo. 

Dans I'exemple de realisation represents sur la figure 2, le filtrage est realise dans le domains frSquentiel. Cette 
implementation est bien adaptee aux tongues reponses impulsionnelles des trajets d'6cho et au degr6 eieve de cor- 

40 relation que presentent les signaux de parole. Les convolutions sont remplacdes par dss produits tsrms k tsrme entrs 
Iss transform6ss de Fourier du signal de reception st ds la reponse impulsionnells, st calcuiess sslon uns tschniqus 
ds convolution rapide OLS (Overiap-and-Save) ou OLA (Overlap-Add) (voir J.J. Shynk : "Frequency-Domain and Mul- 
tirate Adaptive Filtering", IEEE SP Magazine, janvier 1 992, pages 1 5-37). Pour calculsr dss convolutions par un filtrs 
RIF de N coefficients, la technique OLS ou OLA travaille sur des blocs sucessifs de 2N echantillons des signaux, qui 

45 presentent des recouvrements de taille N. 

Les moyens de controle d'adaptation peuvent etre agences pour calculer le gradient selon un algorithme non 
contraint, tel que celui decrit par D. MAN SOUR et al dans "Unconstrained frequency-domain adaptive filter", IEEE 
Trans, on Acoustics, Speech and Signal Processing, Vol. 30, n° 5, Octobre 1982, pages 726-734. Tel est le cas de 
I'exemple represente sur la figure 2. On peut encore faire appel k un algorithme du gradient contraint (voir par exemple 

bo e.R. FERRARA : "Fast implementation of LMS adaptive filters". IEEE Trans, on Acoustics, Speech and Signal Pro- 
cessing, Vol. 28, n° 4, Aout 1980, pages 474-475), les moyens de controle d'adaptation comportant alors des moyens 
de projection pour tenir compte du fait que les 2N coefficients de filtrage ne represented en fait que N termes de la 
reponse impulsionnelle. 

Pour r6duire le retard inherent au traitement par blocs et am6liorer les performances, le filtrage dans le domaine 
55 frequentiel peut etre effectue avec un partitionnement de la reponse impulsionnelle (voir J.S. Soo et al : "Multidelay 
Block Frequency Domain' Adaptive Filtsr", IEEE Trans, on Acoustics, Speech and Signal Processing, Vol.38, N°2, 
f6vrier 1990, pages 373-376). Une autrs possibility est d'employer une technique de recouvrement telle que celle 
d6crite parO. Ait Amrane et al : "Low-Delay Frequency Domain LMS Algorithm", IEEE Proc. ICASSP'92, Vol.4, Mars 
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1992, pages 9-12. 

D'autre part, les annuleurs d'Echo dans le domaine frEquentiel ont I'avantage de bien se prEter.Si des traitements 
da reduction du bruit dans le signal d'erreur (cf. J. BOUDY et al : "A globally optimised frequency<Jomain acoustic 
echo canceller for adverse environment applications", Proc. of the International Workshop on Acoustic Echo and Noise 

s Control, Roros, Norvege, Juin 1995, pages 95-98 ; ou la demande internationale de brevet W096/23384). 

Dans I'annuleur d'Echo de la figure 2, des convertisseurs sErie-parallEle 11, 12 recoivent respectivement les si- 
gnaux x(n) et z(n), et les mettent sous forme de segments successifs non-recouvrants de S Echantillons, S Etant un 
entier infErieur ou Egal a la longueur N de la rEponse impulsionnelle. Le signal de prediction d'Echo y(n) est egalement 
produit sous forme de segments correspondants de S Echantillons. Le soustracteur 8 realise alors foperation e(n)=z 

10 (n)-y(n) en parallEle pour produire le signal d'erreur sous forme de segments de S Echantillons qui sont ensuite remis 
sous forme sEquentielle par un convertisseur parallel e-s6rie 1 3. 

Du fait de la technique de convolution employee, qui est une technique OLS dans le cas represents, les transfor- 
mEes de Fourier des signaux x(n) et e(n) sont rEalisEes sur des blocs de longueur 2N echantillons, par des unites 16, 
17 mettant en oeuvre un algorithme conventionnel de transformation de Fourier rapide (TFR). Les blocs successifs 

is du signal de rEception x(n) sont construits par une unite 1 8 de mani&re k presenter des recouvrements de 2N-S Echan- 
tillons. Ces blocs sont numErotEs par un indice k. A rEception d'un segment de S Echantillons du signal x(n), I'unitE 1 8 
constitue le bloc k en supprimant les S Echantillons les plus anciens du bloc k-1, en dEcalant les 2N-S echantillons 
restants, et en introduisant les S nouveaux Echantillons. Les blocs successifs du signal d'erreur e(n) sont construits 
par une unite 19 k partir des segments dElivrEs par le soustracteur 8 de facon que chaque bloc de longueur 2N com- 

20 mence par N Echantillons k zero et que les N derniers Echantillons de deux blocs successifs prEsentent des recou- 
vrements de N-S echantillons. On note respectivement X^f) et E k (f) (1<f<2N) les composantes de Fourier obtenues 
par transformation des blocs k du signal de rEception x(n) et du signal d'erreur e(n). Pour une f rEquence d'Echantillon- 
nage de 8 kHz, des tallies de segments et de blocs typiques pour le filtrage dans le domaine frEquentiel sont S=32, 
2N=256. 

2$ Les composantes de Fourier produites par les unites 16, 17 sont fournies aux moyens 9 de controle d'adaptation 

qui, dans I'exemple illustre par la figure 2, consistent en 2N modules independants 9.f (1<!<2N) recevant chacun les 
composantes correspondantes X^t) et E k (f) et produisant respectivement les 2N coefficients complexes servant au 
filtrage dans le domaine frEquentiel. Le filtre 7 se compose de 2N multipliers complexes 7.f recevant chacun une 
composante de Fourier X^f) et le coefficient correspondant dElivrE par le module 9.f de contrdle d'adaptation, et dEli- 

30 vrant une composante de Fourier Y k (f). En pratique, compte tenu des propriEtEs de symEtrie de la transtormEe de 
Fourier, on peut obtenir les 2N composantes Y k (f) au moyen de N+1 multiplieurs seulement. Une unite 21 assure une 
transformEe de Fourier rapide inverse (TFRI) du bloc de 2N composantes Y k (f). A partir des blocs de 2N Echantillons 
dElivrEs par I'unitE 21 f une unite 22 construit segment par segment le signal de prEdiction d'Echo y(n). L'unitE 22 extrait 
les N derniers Echantillons de chaque bloc recu de I'unitE 21 et ( lorsque S<N, assure en outre les additions k recou- 
ps vrement nEcessaires a la reconstitution des segments successifs du signal filtre y(n). 

Comme il est bien connu, les transformEes de Fourier rapides peuvent etre remplacEes par des transformEes de 
Hartley rapides, qui conduisent k des composantes spectrales directement liEes aux composantes de Fourier (voir T. 
W. Wong et al : "Adaptive Filtering Using Hartley Transform and Overlap-Save Method", IEEE Trans, on Signal Pro- 
cessing, Vol.39, N 9 7, juillet 1991, pages 1708-1711). 

40 La figure 3 montre une autre forme de realisation possible d'un annuleur d'echo adaptatif, faisant appel k une 

dEcomposition en sous-bandes. On considere k titre d'exemple le cas ou sont prEvues quatre sous-bandes. On note 
x (i) ( n ')» ( n ') ■ z (i) ( n "> et e (D ,ss com P osantes relatives k la sous-bande i ( 1 3£4) des signaux x(n), y(n), z(n) et e(n). 

Les signaux x (i) (n') et z (f) (n') sont obtenus en dEcomposant les signaux x(n) et z(n) au moyen de deux Etages de 
dEcomposition, le premier Etage comprenant un module D, D' k une entrEe et deux sorties du type reprEsentE sur la 

45 figure 4, et le second Etage comprenant deux modules D, D' du m&me type dont les entrEes sont respectivement 
reliEes aux deux sorties du module du premier Etage. Chaque module D ou D* comprend un filtre passe^bas 23 et un 
filtre passe-haut 24 recevant tous deux le signal d'entrEe du module et deux circuits 26 qui sous-Echantillonnent dans 
un facteur 2 les signaux de sortie des filtres 23, 24, pour produire les deux signaux de sortie du module D, D\ 

L'annuleur d'Echo en sous-bandes comporte un filtre adaptatif 7 pour chaque sous-bande, qui filtre dans le domaine 

50 tempore! la composante correspondante x (i) (n') du signal de rEception pour produire la composante y (i) (n') du signal 
de prEdiction d'Echo qui est soustraite par un soustracteur B de la composante z (l) (n 1 ) du signal d'observation. Pour 
chaque sous-bande, un module 9 controle I'adaptation du filtre 7. 

Des moyens de recombinaison en deux Etages construisent le signal d'erreur e(n) k partir des quatre composantes 
e rn (n') dElivrEes par les soustracteurs 8. Le premier Etage comprend deux modules R k deux entrEes et une sortie du 

55 type reprEsentE sur la figure 5, et le second Etage comprend un seul module R du meme type dont les deux entrEes 
sont reliEes aux sorties des deux modules du premier Etage. Chaque module R comprend deux circuits 27 de surE- 
chantillonnage dans un facteur 2 des deux signaux d'entrEe du module, un filtre passe-bas 28 et un filtre passe-haut 
29 recevant respectivement les deux signaux surEchantillonnEs produits par les circuits 27, et un sommateur 31 qui 



EP 0 854 626 A1 



additionne les signaux de sortie des deux filtres 28, 29 pour produire le signal de sortie du module R. 

On connait diverses formes classiques pour les filtres de decomposition et de recombinaison 23, 24, 28, 29 qui 
peuvent gen6ralement etre a reponse impulsionnelle finie ou infinie, par exemple des filtres conjugues en quadrature 
(CQF), des filtres miroirs en quadrature (QMF), ou encore des filtres a decomposition en ondelettes (WDF). D'autre 

5 part, on peut utiliser, plutot que la structure en plusieurs etages, une structure a decomposition directe de type banc 
de filtres modules en cosinus (voir H.S. Matvar : "Extended Lapped Transforms : Properties, Applications, and Fast 
Algorithms - , IEEE Trans, on Signal Processing, Vol.40, N°11, novembre 1992, pages 2703-2714, ou P.P. 
Vaidyanathan : "Multirate Digital Filters, Filter Banks, Polyphase Networks and Applications : A Tutorial", Proc. of the 
IEEE, Vol.78, N°1, janvier 1990, pages 56-93). 

10 Les algorithmes d'adaptation mis en oeuvre par les modules 9 de I'annuleur d'echo en sous-bandes peuvent etre 

du meme type que ceux couramment utilises pour adapter le filtrage tempore! en pleine bande. II existe neanmoins 
des algorithmes particulierement bien adaptes a la structure en sous-bandes, tels que I'algorithme rapide des moindres 
carres recursifs (RLS) a decomposition QR (cf. F. CAPMAN et al : "Acoustic echo cancellation using a fast-QR-RLS 
algorithm and multirate schemes", IEEE ICASSP'95, Detroit, 1995; ou la demande intemationale de brevet 

T5 W096/21313). II est egalement possible d'utiliser des parametres de convergence et/ou des algorithmes d'adaptation 
difterents d'une sous-bande a une autre. On peut notamment utiliser dans les sous-bandes relatives aux frequences 
les plus basses un algorithme (tel que le QR-HLS rapide) plus emcace que dans les sous-Danaes supeneures, ou on 
peut utiliser un algorithme moins exigeant en termes de ressources de calcul (tel qu'un NLMS). 

Dans certains annuleurs d'echo, les coefficients du filtrage de prediction d'echo sont figes lors du fonctionnement. 

20 Ces coefficients peuvent etre fixes ; ils sont alors determines une fois pour toutes de maniere a prendre en compte, 
par exemple, les reponses du haut-parleur 3 et du microphone 4, les couplages directs entre le haut-parleur et le 
microphone, et eventuellement les couplages indirects si le haut-parleur et le microphone restent toujours dans les 
memes positions. Les coefficients peuvent encore §tre determines par apprentissage, soit a la mise en marche de 
I'annuleur d'echo soit periodiquement, en faisant converger un algorithme d'adaptation sur un signal superpose ou 

25 substitue au signal de reception et ayant des proprietes de decorrelation analogues a celles d'un bruit blanc. 

On connait egalement des annuleurs d'echo dans lesquels il y a plusieurs filtres de prediction d'echo, le signal 
emis vers le canal de communication 6tant obtenu en s6lectionnant I'un des filtres selon des criteres predefinis tels 
qu'une minimisation de Penergie du signal d'erreur, ou encore en combinant les sorties de plusieurs filtres de prediction 
d'echo (voir W. Armbruster : "Wideband Acoustic Echo Canceller with Two Filter Structure", Signal Processing VI : 

30 Theories and Applications, 1992, pages 1611-1614, ou J.S. Soo et al : "A Multistep Size (MSS) Frequency Domain 
Adaptive Filter", IEEE Trans, on Signal Processing, Vol.39, n°1, janvier 1991, pages 115-121). L'annuleur d'echo peut 
par exemple comporter un ou plusieurs filtres adaptatifs fonctionnant avec des parametres de convergence diff6rents 
ou selon des algorithmes d'adaptation differents, et/ou un ou plusieurs filtres a coefficients fixes, et/ou un ou plusieurs 
filtres dont les coefficients sont determines par apprentissage. 

3S Un but principal de la presente invention est d'ameiiorer la convergence et/ou la qualite du signal de sortie d'un 

annuleur d'echo, notamment d'un annuleur d'echo faisant intervenir une certaine resolution spectrale (filtrage dans le 
domaine f requentiel ou en sous-bandes). 

L'invention propose ainsi un precede d'annulation d'echo, dans lequel on applique un signal de reception en entr6e 
d'un systeme generateur d'echo, on recueille un signal d'observation en sortie dudit systeme, on filtre le signal de 

40 reception pour produire au moins un signal de prediction d'un echo du signal de reception dans le signal d'observation, 
et on soustrait un signal de prediction d'echo du signal d'observation pour produire un signal d'erreur. Selon l'invention, 
on decompose le signal d'observation en M composantes spectrales se rapportant a des frequences differentes du 
spectre, M designant un nombre entier plus grand que 1 , on obtient M composantes spectrales correspondantes d'un 
signal de prediction d'6cho, on calcule M grandeurs a partir des composantes spectrales du signal d'observation et du 

45 signal de prediction d'6cho, chacune desdites M grandeurs ayant une valeur a variation progressive representative 
d'un degr6 de similitude entre, respectivement, I'une des M composantes spectrales du signal d'observation et la 
composante spectrale correspondante du signal de prediction d'echo, et on opere au moins une mise en forme spec- 
trale du signal d'erreur en fonction desdites M grandeurs calcuiees. 

Les M grandeurs permettent de disposer d'information relative a la presence d'interferences (bruit ambiant ou 

so parole locale) dans le signal d'observation, cette information etant selective en frequence et a variation progressive, 
ce qui assure une richesse d'information et une precision sensiblement plus importantes que les detecteurs de double 
parole classiques qui sont essentiellement a decision binaire et sur toute la bande spectrale. 

Lorsqu'on adapte le filtrage du signal de reception en fonction du signal de reception et du signal d'erreur, la mise 
en forme spectrale du signal d'erreur peut servir a produire un signal d'erreur filtre utilise pour I'adaptation du filtrage 

55 du signal de reception. La mise en forme consiste alors a attenuer, de preference de tacon non lineaire, les compo- 
santes spectrales du signal d'erreur pour lesquelles ta grandeur calcul6e traduit une similitude relativement faible entre 
les composantes spectrales correspondantes du signal d'observation et du signal de prediction d'6cho. 

La mise en forme spectrale du signal d'erreur peut egalement servir a produire un signal d'erreur post-filtr6 que 
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I'annuleur d'echo 6met vers le canal de communication. Dans ce cas, le post-filtrage tend plutet a attenuer, de pr6f6- 
rence de facon non Iin6aire, les composantes du signal d'erreur pour lesquelles la grandeur calcuiee montre une 
similitude relativement importante entre les composantes correspondantes des signaux y(n) et z(n). 

Contrairement a la pratique habituelle dans le domaine des detecteurs de double parole, ce sont les composantes 

5 spectrales du signal de prediction d'echo et non celles du signal de reception qui sont compares aux composantes 
correspondantes du signal d'observation pour quantifier le degre de similitude entre ces composantes. Ainsi, les as- 
pects temporels des trajets d'echo sont deja pris en compte par le filtrage op6r6 sur le signal de reception. 

On note que si I'annuleur d'echo est pourvu de plusieurs filtres de prediction d'echo, il est concevable que le signal 
de prediction d'echo sur la base duquel sont calculees les M grandeurs representant les degr6s de similitude ne soit 

10 pas le meme que celui qui est soustrait du signal d'observation pour produire le signal d'erreur soumis a la mise en 
forme spectrale (filtrage et/ou post-filtrage). On peut notamment prevoir qu'un filtre a coefficients fixes soit employe 
sp6cialement pour produire le signal de prediction d'echo servant au calcul des M grandeurs de similitude. 

Lorsque les M composantes spectrales du signal d'observation et du signal de prediction d'6cho sont des compo- 
santes de Fourier de ces signaux, il est avantageux que la grandeur representative du degr6 de similitude entre deux 

is composantes de Fourier correspondantes du signal d'observation et du signal de prediction d'6cho soit une valeur 
d'unefonction de coherence entre ces composantes de Fourier. Cette disposition pourraavantageusement 6tre retenue 
dans un annuleur d'6cho a filtrage dans le domaine fr6quentiel, puisque des calculs de transform6es de Fourier ou de 
Hartley sont de toute facon necessaires dans ce cas. Ces dispositions sont 6galement applicables dans le cas d'un 
filtrage dans le domaine temporel. 

20 L'utilisation de la fonction de coherence a d6ja ete proposee dans des detecteurs de double parole (cf . T. G ANSLER 

et al : "A double-talk detector based on coherence", Proc. of the International Conference on Signal Processing Ap- 
plications and Technology, Boston, Octobre 1995, pages 332-336). Outre le contexte duplication different, on peut 
noter que ces auteurs calculent la coherence entre le signal d'observation et le signal de reception, et done ne corrigent 
pas le canal acoustique, et que les decisions prises en relation avec les coherences calculees sont de nature binaire. 

2$ En variante, il est possible de remplacer la fonction de coherence par une fonction decroissante de la valeur 

absolue de I'angle lorme, dans le plan complexe, entre les composantes de Fourier correspondantes du signal d'ob- 
servation et du signal de prediction d'echo. Cette fonction decroissante peut notamment etre le cosinus, ou le cosinus 
carr6, de I'angle bom6 a±n/2. On peut montrer que ce cosinus est directement lie a la fonction de coherence. L'avantage 
de cette formulation est que le cosinus, ou son carre, peut 6tre calcuie simplement sans avoir besoin d'une m6moire 

30 necessaire aux calculs de moyennes utilises pour obtenir les valeurs de la fonction de coherence. La mesure reflete 
ainsi plus rapidement les variations des signaux. 

Dans le cas d'un annuleur d'echo en sous-bandes, les M composantes spectrales du signal d'observation sont 
obtenues par une decomposition en M sous-bandes du signal d'observation, la meme decomposition en sous-bandes 
est effectu6e sur le signal de reception pour produire M composantes spectrales correspondantes du signal de r6cep- 

35 tion, le signal de prediction d'echo est produit sous la forme de ses M composantes spectrales respectivement deiivrees 
par M filtres adaptatifs auxquels sont respectivement appliqu6es les M composantes spectrales du signal de reception, 
et le signal d'erreur est produit sous la forme de M composantes spectrales correspondantes dont chacune est obtenue 
en soustrayant une composante spectrale du signal de prediction d'echo d'une composante spectrale correspondante 
du signal d'observation. Dans ce cas, la grandeur representative du degre de similitude entre deux composantes 

40 spectrales correspondantes du signal d'observation et du signal de prediction d'6cho est avantageusement une valeur 
d'une fonction d'intercorr6lation normalis6e entre ces composantes spectrales. 

D'autres particularites et avantages de la pr6sente invention apparaitront dans la description ci-apres d'exemples 
de realisation non limitatifs, en reference aux dessins annexes, dans lesquels : 

45 . la figure 1, pr6c6demment comment6e, est un schema de principe d'un annuleur d'6cho adaptatif a filtrage en 
pleine bande dans le domaine temporel ; 

- la figure 2, prec6demment commentee, est un schema de principe d'un annuleur d'6cho adaptatif a filtrage dans 
le domaine f requentiel ; 

- la figure 3, precedemment commentee, est un schema de principe d'un annuleur d'6cho adaptatif a filtrage en 
so sous-bandes ; 

- la figure 4 et 5, pr6cedemment commentees, sont des schemas synoptiques de moyens de decomposition et de 
recombinaison de I'annuleur d'echo de la figure 3 ; 

la figure 6 est un schema de principe d'un annuleur d'echo adaptatif selon la presente invention ; 

- la figure 7 est un schema synoptique d'une unite de calcul de coherence de I'annuleur d'6cho de la figure 6 ; 

55 . la figure 8 est un schema synoptique de moyens de filtrage du signal d'erreur dans un annuleur d'echo selon la 
figure 6 ; 

- les figures 9 et 10 sont des graphiques illustrant des formes possibles de fonctions non Iin6aires utilisables pour 
mettre en oeuvre I'invention ; 
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les figures 11 et 1 2 sont des sch6mas de variantes de realisation des moyens de filtrage de la figure 8 ; et 

la figure 1 3 est un schema synoptique de moyens permettant la mise en oeuvre de la pr6sente invention dans un 

annuleur d'6cho adaptatif k filtrage en sous-bandes. 

5 La figure 6 reprend le schema de principe de I'annuleur cfecho de la figure 1 , en y ajoutant trois modules 32, 33, 

34 servant k la mise en oeuvre de la prSsente invention. Le module 32 calcule M grandeurs qui, avec une certaine 
resolution spectrale, caracterisent un degre de similitude entre le signal d'observation z(n) et le signal de prediction 
d'echo y(n). Le module 33 est un filtre qui, sur la base de ces M grandeurs, opere une mise en forme spectrale du 
signal d'erreur e(n), et fournit le signal ainsi mis en forme e1(n) au module de controle d'adaptation 9. Le module 34 

10 est un post-filtre qui opere une autre mise en forme spectrale du signal d'erreur e(n) sur la base des M grandeurs pour 
emettre le signal ainsi post-filtre* e2(n) vers le canal de communication auquel est relie I'annuleur d'echo. 

Bien que le schema de la figure 6 ait ete etabli dans le cas ou, selon la figure 1 , 1'annulation d'6cho repose sur un 
filtrage en pleine bande dans le domaine temporel, on comprendra que I'adjonction de trois modules tels que 32, 33, 
34 permet de facon g6n6rale de mettre en oeuvre ('invention dans tout type d'annuleur d'echo. 

is Dans le cas de la figure 6, les M grandeurs calcul6es par le module 32 sont des valeurs d'une fonction de coherence 

entre des composantes de Fourier du signal d'observation z(n) et des composantes de Fourier correspondantes du 

signal de prediction d'echo y(n). Le module 32 peut alors avoir la constitution representee sur la figure 7. 

Dans le module 32 de la figure 7, des convertisseurs serie-parall&le 36, 37 regoivent respectivement les signaux 
z(n) et y(n), et les mettent sous forme de blocs successifs de M 6chantillons. Ces blocs peuvent gen6ralement etre 

20 recouvrants ou non-recouvrants. Les transformees de Fourier des signaux y(n) et z(n) sont realis6es sur ces blocs de 
longueur M echantillons par des unites 38, 39 mettant en oeuvre un algorithrne conventionnel de TFR. On note Zp(m) 
et Yp(m) les m-iSmes composantes de Fourier deiivrees par les unites 38 et 39 relativement au bloc de rang p du 
signal d'observation z(n) et du signal de prediction d'echo y(n) (1<m^M). 

Pour chaque index frequentiel m, les composantes de Fourier Zp(m) et Y p (m) sont fournies k une unite associee 

2S 41 .m qui se charge de calculer la valeur correspondante de la fonction de coherence msCp(m). Ces unites 41 .m ont 
toutes la meme constitution. Dans I'exemple considere, la fonction de coherence employee est le module carre de la 
coherence ("magnitude squared- coherence"), qui s'exprime par : 

msc (m) = — ^ — ( 1 ) 

ou (m) et ^(m) d6signent les densit6s spectrales de puissance, estimees au moment du bloc p pour la m-ieme 
3S frequence, des signaux y(n) et z(n), et Y p z (™) designe la densite spectrale croisee de ces deux signaux estimee au 
moment du bloc p pour la m-ieme frequence. Ces densites spectrales de puissance peuvent etre calcul6es simplement 
par le module 41 .m h I'aide d'un facteur d'oubli X tel que 0«X<1 , de la maniere suivante : 



40 



45 



y ™{m) = \.y%{mW-\)\ Y p (m)\ 2 (2) 
yf(m) = Vyf t (mW1-X).| 2 [mf O) 



y Y p z (m) = X.y£(m) + 0-X). Y p (m)[ z p (m)\* (4) 



ou I.I 2 designe le module carre et (.)* le conjugue d'un nombre complexe. 

Pour chaque indice de composante spectrale m, la grandeur msc p (m) presente une variation progressive, croissant 
50 entre 0 el 1 , en fonction du degre de similitude entre les composantes de Fourier Zp(m) et Y p (m). Cette grandeur est 
mise k jour k chaque nouveau bloc de M echantillons. 

La figure 8 montre un schema synoptique du filtre 33 et du post-filtre 34 qui traitent le signal d'erreur e(n). Ces 
traitements sont ici effectues dans le domaine frequentiel, les filtres 33 et 34 partageant les circuits produisant les 
composantes de Fourier E p (m) du signal d'erreur e(n). Ces circuits comprennent un convertisseur s6rie-parallele 53 
55 produisant des blocs de M echantillons du signal d'erreur e(n) de la m6me maniere que les convertisseurs 36, 37 de 
la figure 7, et une unite de calcul de TFR 54.' 

Le filtre 33 comporte M circuits 56.m (1£m£M) calculant chacun, pour une valeur respective msc p (m) de la fonction 
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de coherence, la valeur prise par une fonction non lineaire croissants FNL1. La fonction non lineaire FNL1 a par 
exemple la forme representee sur la figure 9, definie par trois parametres a1,b1 etT1 telsque0<a1<b1<1 etO<T1<1 : 
la fonction FNL1(x) vaut T1 pour CKx£a1, elle croTt de T1 a 1 quant x croTt de a1 ab1, etelle vaut 1 pourb1<x<1. Le 
jeu de parametres T1 , a1 , b1 peut §tre identique pour chacune des M composantes de Fourier, ou il peut etre different 

5 d'une composante de Fourier a une autre. Les valeurs retenues pour ces parametres peuvent etre determinees ex- 
perimentalement ou par simulation en fonction de ('application particuliere et des caracteristiques du terminal equipe 
de I'annuleur d'echo. Chaque valeur FNL1(msc p (m)) produite par un circuit 56.m est foumie a un multiplieur respectif 
57.m recevant par ailleurs la composante de Fourier correspondante E p (m) du signal d'erreur. Les valeurs E1 p (m) 
delivrees par les M multiplieurs 57.m sont fournies a une unite 58 de calcul de TFRI qui obtient leur transformee de 

10 Fourier inverse sous la forme d'un bloc de M echantillons. Un convertisseur parallele-serie 59 forme le signal d'erreur 
filtre el(n) sur la base des blocs delivres par I'unite 58. 

Ce signal e1 (n) est exploite par le module de controle d'adaptation 9 de la meme maniere que le signal e(n) dans 
le cas classique de la figure 1. L'adaptation du filtre 7 peut etre realise selon tout algorithme approprte. 

Le post-filtre 34 comporte M circuits 61. m (1<m<M) calculant chacun, pour une valeur respective msc p (m) de la 

is fonction de coherence, la valeur prise par une fonction non lineaire decroissante FNL2. La fonction non lineaire FNL2 
a par exemple la forme representee sur la figure 10, definie par trois parametres a2, b2 et 12 tels que 0<a2<b2<1 et 
0<T2<1 : la fonction FNL2(x) vaut 1 pour C£x£a2, elle decroTt de 1 a T2 quant x croTt de az a 02, et elle vaui T2 pour 
b2<x^1. Le jeu de parametres T2, a2, b2 peut etre identique pour chacune des M composantes de Fourier, ou il peut 
etre different d'une composante de Fourier a une autre. Les valeurs retenues pour ces parametres peuvent etre de- 

20 termin<§es experimentalement ou par simulation. Chaque valeur FNL2(msc p (m)) produite par un circuit 61 .m est foumie 
a un multiplieur respectif 62.m recevant par ailleurs la composante de Fourier correspondante E p (m) du signal d'erreur. 
Les valeurs E2 p (m) delivrees par les M multiplieurs 62.m sont fournies a une unite 63 de calcul de TFRI qui obtient 
leur transformee de Fourier inverse sous la forme d'un bloc de M echantillons. Un convertisseur parallele-serie 64 
forme le signal d'erreur post-filtre e2(n) sur la base des blocs delivres par I'unite 63. 

2S Lorsqu'il est prevu des modules tels que 32, 33, 34 (figure 6) dans le cadre d'un annuleur d'echo fonctionnant par 

filtrage dans le domaine frequentiel, il est a noter que la resolution spectrale utilises pour les calculs de coherence 
peut etre egale ou inferieure a la resolution spectrale du filtre de moderation des echos. En d'autres termes, avec 
les notations de la figure 2, on peut avoir M=2N ou M<2N. D'autre part, la cadence de mise a jour des calculs de 
coherence, qui correspond a la cadence de renouvellement des blocs p de M echantillons sur la figure 7, peut etre 

30 differente de la cadence de renouvellement des blocs k de 2N echantillons sur la figure 2. 

En outre, il est possible de subdiviser le spectre des signaux en Q sous-bandes (Q>1) pour les filtrages du signal 
d'erreur. Les valeurs de la fonction de coherence pour les composantes frequentielles relevant d'une sous-bande 
donnee sont alors moyennees, et c'est la valeur moyenne ainsi obtenue qui est soumise aux fonctions non-lineaires 
FNL1 et FNL2. La valeur d'une fonction non-lineaire pour une sous-bande donnee sert alors a ponderer chacune des 

35 composantes de Fourier du signal d'erreur qui relevent de cette sous-bande. 

Les figures 11 et 12 montrent de telles implementations pour le filtre et le post-filtre du signal d'erreur dans le cas 
d'un annuleur d'echo travaillant dans le domaine frequentiel conlormement a la figure 2, les valeurs de la fonction de 
coherence msc p (m) (1<m<M avec M<2N) etant obtenues par un module 32 conforme a ceiui represents sur la figure 
7. On suppose que la q-ieme sous-bande (1<Sq£Q) couvre M q composantes de Fourier des signaux z(n) et y(n) pour 

40 lesquelles le module 32 a calcule les valeurs de la fonction de coherence, et couvre N q composantes de Fourier du 
signal d'erreur calculees par I'unite de TFR 17. La subdivision en Q sous-bandes peut etre unifomne (M 1 =M 2 =... =M Q 
et N 1= N 2 =... =N Q ) ou non uniforme, par exemple selon une echelle de Bark telle qu'utilisee dans le domaine de la 
psycho-acoustique. 

Le filtre du signal d'erreur represents sur la figure 11 comprend Q circuits d'integration 66 q (1£q£Q) qui calculent 
4$ chacun la moyenne arithm6tique des M q valeurs de la fonction de coherence relatives a la sous-bande q, et Q circuits 
67.q qui calculent, pour les moyennes delivrees par les circuits d'integration 66q, les valeurs prises par la fonction 
non-lineaire croissante FNL1 qui peut etre identique ou differente pour chaque sous-bande. Le filtre comprend en outre 
2N multiplieurs 68.f (1<f<2N) recevant chacun une composante de Fourier respective E k (f) du signal d'erreur ainsi que 
la valeur delivree par le circuit 67q, ou q designe la sous-bande dont releve la composante de Fourier f du signal 
50 d'erreur. Dans le cas de la figure 11, le signal d'erreur filtre e1(n) n'est pas produit explicitement dans le domaine 
temporel, mais sous la forme de ses 2N composantes de Fourier E 1 k (f) qui sont les composantes de Fourier ponderees 
du signal d'erreur respect ivement delivrees par les multiplieurs 68.f. Ces composantes E1 k (f) sont respectivement 
fournies en tant que signal d'erreur filtre aux 2N modules de controle d'adaptation 9.f qui sont du meme type que ceux 
represents sur la figure 2. 

55 Les Q circuits d'integration 66.q font egalement partie du post-filtre du signal d'erreur represente sur la figure 12, 

qui comprend Sgalement Q circuits 69.q qui calculent, pour les moyennes delivrees par les circuits d'integration 66.q, 
les valeurs prises par la fonction non-lineaire decroissante FNL2 qui peut etre identique ou differente pour chaque 
sous-bande. Le post-filtre comprend en outre 2N multiplieurs 7 1 .f (1£f<£N) recevant chacun une composante de Fourier 
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respective E k (f) du signal d'erreur ainsi que la valeur d6livred par le circuit 69.q, ou q dSsigne la sous-bande dont 
reieve la composante de Fourier f du signal d'erreur. Les composantes de Fourier ponderers E2 k (f) respectivement 
delivreds par les multipliers 68.f sont f ournies a une unite 72 de calcul de TFRI produisant des blocs de 2N echantillons. 
A partir de ces blocs, une units 73 semblable a I'unite 22 de la figure 2 construit segment par segment le signal d'erreur 

5 post-filtr§ e2(n) qu'un convertisseur parallele-s§rie 74 remet sous forme sequentielle. 

Comme il a et§ indiquS plus haul, une variante avantageuse de ^invention consiste a remplacer la fonction de 
coherence entre les composantes Zp(m) et Y p (m) par une fonction decroissante de Tangle <p p (m) forme\ dans le plan 
complexe, entre les nombres complexes Zp(m) et Y p (m) qui, dans le cas illustre" par la figure 7 sont produits par les 
unites de TFR 38 et 39. En d'autres termes, Tangle cp p (m) est Targument du nombre complexe Zp(m)/Y p (m). Si on note 

10 a et a' les parties replies respectives de Zp(m) et Y p (m), et b et b' leurs parties imaginaires, alors : 

/ x aa'+bb' 
coscp (m)= == 

J(a %b 2 )(a%b*) 

is 

Ce cosinus de Tangle <p p (m), born6 entre -n/2 et W2 (car un d6phasage de n denote une predence d'intendrences 

etant donne qu'une compensation du canal acoustique est d6ja opered dans le signal y(n) par le nitre (), esi une auire 
mesure normalised du degr6 de similitude entre les composantes Zp(m) et Y p (m), qui predente Tavantage de pouvoir 
etre calculee aisement bloc par bloc sans avoir a memoriser des quantites pour redliser des moyennages (selon les 

20 relations (2) a (4)). En pratique, on pourra calculer le carre ctp(m) de cos<p p (m) pour evrter le calcul de racine carree 
qui pourra etre pris en compte dans la definition des fonctions non Iin6aires FNL1 et FNL2. 

Dans cette variante, les unites 41 .m de la figure 7 sont agencees pour delivrer cc p (m)=0 si aa'+ bb'^), et a p (m) = 
(aa'+ bb') 2 /[(a 2 + b 2 )(a' 2 + b' 2 )] sinon. Ces grandeurs cc p (m) peuvent etre traitees de facon semblable aux msc p (m) par 
des filtres et post-filtres tels que ceux schematises sur les figures 9 a 12. 

2S Dans le cas d'un annuleur d'echo en sous-bandes, la resolution spectrale dont on a besoin selon Invention ne 

resulte pas d'une transformed de Fourier, mais d'une decomposition en M sous-bandes. Les M grandeurs servant a 
quantifier de maniere progressive la similitude entre les signaux z(n) et y(n) peuvent alors etre les valeurs prises par 
une fonction d'intercorreMation normalised entre les composantes spectrales (n') et y (i) (n*) (figure 3), telle que : 
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ou E[.] designe Tesp6rance mathematique. Les esperances mathematiques de la relation (5) peuvent etre calculees 
par moyennage exact sur une fenStre d'analyse, auquel cas les r (j) sont mis a jour a chaque nouvelle fenStre, ou au 
40 moyen d'un facteur d'oubli de facon semblable aux relations (2) a (4), auquel cas les r 0) peuvent etre mis a jour a 
chaque temps de sous-echantillonnage. La valeur de r (i) croTt progressivement de 0 a 1 selon le degre de similitude 
entre les signaux z (l) (n') et y (i) (n'). 

Chaque voie de traitement d'une sous-bande i peut alors etre completed par les moyens represents sur la figure 

13, qui comprennent : 

45 

- une unit6 76 calculant la valeur r (j) de la fonction d'intercor relation normalised pour les composantes z (i) (n') et y (j) 
(n') fournies respectivement par les moyens de decomposition D' et par le filtre 7 relatif a la sous-bande i ; 

- un circuit 77 calculant pour la valeur de r (i) la valeur prise par une fonction non lineaire croissante FNL1 qui peut 
etre du meme type que celle represented graphiquement sur la figure 9 et qui peut etre difterente d'une sous- 

so bande a une autre ; 

- un multiplier 78 qui pondere la composante spectrale e^ (n l ) du signal d'erreur delivree par le soustracteur 8 par 
la valeur foumie par le circuit 77, et qui fournit la composante pondered e1 (i) (n') au module de controle d'adaptation 
9 relatif a la sous-bande i ; 

- un circuit 79 calculant pour la valeur de r w la valeur prise par une fonction non linedire decroissante FNL2 qui peut 
55 etre du meme type que celle representee graphiquement sur la figure 10 et qui peut etre differente d'une sous- 
bande a une autre ; et 

- un multiplier 81 qui pondere la composante spectrale e^n') du signal d'erreur delvree par le soustracteur 8 par 
la valeur foumie par le circuit 79, et qui foumit la composante pond6ree e2 (0 (n') au module R approprte des moyens 
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de recombinaison. 

Les M circuits 77 et les M multiplieurs 78 relatifs aux diff6rentes sous-bandes constituent, de facon semblable au 
filtre 33 de la figure 6, le filtre qui opere ta mise en forme spectrale du signal d'erreur en fonction des M grandeurs r (i) 
5 pour ('adaptation du filtrage du signal de reception x(n). 

Les M circuits 79, les M multiplieurs 81 et les modules R de recombinaison constituent, de facon semblable au 
post-filtre 34 de la figure 6, le post-filtre qui opere la mise en forme spectrale du signal d'erreur en fonction des M 
grandeurs pour la production du signal de sortie de I'annuleur d'echo. 

Pour le calcul des esperances math<§matiques de la relation (5), I'unite 76 peut comprendre trois circuits de moyen- 
10 nage 82, 83, 84 recevant respectivement les valeurs successives de y (i) (n').z (i ) (n'), de y 2 (i) (n') et de z 2 (i) (n 1 ) d'un 
multiplieur 86 et de deux circuits de mise au carre 87, 88. Le quotient de la relation (5) est calcule par un diviseur 89 
recevant le numSrateur d'un circuit de mise au carre 91 et le denominateur d'un multiplieur 92. 

Dans la description ci-dessus et les dessins annexes, on a presente des modes de realisation de ('invention a 
partir d'agencements matdriels de modules, d'unites ou de circuits individualists, afin de rendre I'explication aussi 
is claire que possible. L'homme du metier appr6ciera que I'invention peut 6galement etre mise en oeuvre diff6remment, 
notamment en programmant un processeur de signal num6rique approprie ainsi qu'il est usuel dans la technique des 
annuleurs d'echo. 



20 Revendications 

1. Procede d'annulation d'echo, dans lequel on applique un signal de reception (x(n)) en entree d'un systeme gene- 
rates d'6cho (1 ), on recueille un signal d'observation (z(n)) en sortie dudit systeme, on filtre le signal de reception 
pour produire au moins un signal (y(n)) de prediction d'un echo du signal de reception dans le signal d'observation, 

25 et on soustrait un signal de prediction d'echo du signal d'observation pour produire un signal d'erreur (e(n)), ca- 

racterise en ce qu'on decompose le signal ^observation (z(n)) en M composantes spectrales (Zp(m) ; z^ (n')) se 
rapportant a des frequences differentes du spectre, M designant un nombre entier plus grand que 1 , on obtient M 
composantes spectrales correspondantes (Y p (m) ; y 0) (n')) d'un signal de prediction d'echo (y(n)), on calcule M 
grandeurs (msc p (m) ; a p (m) ; r (i) ) a partir des composantes spectrales du signal d'observation et du signal de 

30 prediction d'echo, chacune desdites M grandeurs ayant une valeur a variation progressive representative d'un 

degr6 de similitude entre, respectivement, I'une des M composantes spectrales du signal d'observation et la com- 
. posante spectrale correspondante du signal de prediction d'6cho, et on opere au moins une mise en forme spectrale 
d'un signal d'erreur (e(n)) en fonction desdites M grandeurs calcul6es. 

35 2. Procede selon la revendication 1 , dans lequel on on adapte le filtrage du signal de reception en fonction du signal 
de reception et d'un signal d'erreur, et dans lequel on opere une mise en forme spectrale du signal d'erreur (e(n)) 
produisant un signal d'erreur filtre (e1(n)) en fonction duquel est effectuee I'adaptation du filtrage du signal de 
reception (x(n)). 

40 3. Precede selon la revendication 1 ou 2, dans lequel on opere une mise en forme spectrale d'un signal d'erreur (e 
(n)) produisant un signal d'erreur post-filtre (e2(n)) emis vers un canal de communication. 

4. Procede selon I'une quelconque des revendications 1 a 3, dans lequel les M composantes spectrales (Zp(m),Y p 
(m)) du signal d'observation (z(n)) et du signal de prediction d'6cho (y(n)) sont des composantes de Fourier de 

45 ces signaux, et dans lequel la grandeur (msc p (m)) representative du degr6 de similitude entre deux composantes 

de Fourier correspondantes (Zp(m), Y p (m)) du signal d'observation et du signal de prediction d'6cho est une valeur 
d'une fonction de coherence entre ces composantes de Fourier. 

5. Procede selon I'une quelconque des revendications 1 a 3, dans lequel les M composantes spectrales (Zp(m),Y p 
50 (m)) du signal d'observation (z(n)) et du signal de prediction d'echo (y(n)) sont des composantes de Fourier de 

ces signaux, et dans lequel la grandeur (Op(m)) representative du degr6 de similitude entre deux composantes de 
Fourier correspondantes (Zp(m) ,Y p (m)) du signal d'observation et du signal de prediction d'echo est une valeur 
d'une fonction decroissante de la valeur absolue de Tangle (<p p (m)) forme, dans le plan complexe, entre ces com- 
posantes de Fourier. 

55 

6. Procede selon la revendication 4 ou 5, dans lequel on obtient M composantes de Fourier (E p (m)) du signal d'erreur 
(e(n)) qui correspondent respectivement aux M composantes de Fourier (Zpfm^Ypfm)) du signal d'observation (z 
(n)) et du signal de prediction d'6cho (y(n)), et dans lequel chaque mise en forme spectrale du signal d'erreur (e 
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(n)) comporte une ponderation de chacune de ses M composantes de Fourier (E p (m)) par une quantite ayant une 
dependance non lineaire de la valeur de celfe des M grandeurs (msCp(m) ; Op(m)) qui represente le degre de 
similitude entre les composantes de Fourier correspondantes (Zp(m),Y p (m)) du signal d'observation et du signal 
de prediction d'6cho. 

5 

7. Procede selon la revendication 4 ou 5, dans lequel on obtient un nombre 2N au moins egal a M de composantes 
de Fourier (E k (f)) du signal d'erreur (e(n)), dans lequel on repartit dans plusieurs sous-bandes predeterminees les 
composantes de Fourier des signaux, et on calcule pour chaque sous-bande une moyenne de celles des M gran- 
deurs (mscp(m) ; Op(m)) qui represented des degres de similitude entre les composantes de Fourier relatives a 
10 ladite sous-bande du signal d'observation et du signal de prediction d'echo, et dans lequel chaque mise en forme 

spectrale du signal d'erreur (e(n)) comporte une ponderation de chacune de ses 2N composantes de Fourier (E|< 
(f)) par une quantite ayant une dependance non lineaire de la moyenne calculee pour la sous-bande dont releve 
ladite composante de Fourier du signal d'erreur. 

is 8. Procede selon la revendication 2 et Tune quelconque des revendications 6 et 7, dans lequel, pour la production 
du signal d'erreur filtre (e1(n)), ladite dependance non lineaire est definie par une fonction non lineaire croissante 
(FNL1). — " " 

9. Procede selon la revendication 8, dans lequel on produit le signal d'erreur filtre (e1(n)) en appliquant une trans- 
20 formee de Fourier inverse aux composantes de Fourier ponderees (E1 p (m)) du signal d'erreur. 

10. Procede selon la revendication 8, dans lequel on effectue une transformee de Fourier du signal de r6ception (x 
(n)) pour en obtenir un nombre 2N, au moins egal a M, de composantes de Fourier (X^f)), et on produit un signal 
de prediction d'echo (y(n)) sous forme de 2N composantes de Fourier (Y k (f)) obtenues en ponderant les compo- 

25 santes de Fourier (X k (f)) du signal de reception par des coefficients complexes adaptatifs respectifs, chacun de 

ces coefficients etant associS a une composante de Fourier ponderee (E 1 k (f )) du signal d'erreur qui sert, avec la 
composante de Fourier correspondante (X^f)) du signal de reception, a adapter la valeur dudit coefficient. 

Procede selon la revendication 3 et I'une quelconque des revendications 6 et 7, dans lequel, pour la production 
du signal d'erreur post-filtre (e2(n)), ladite dependance non lineaire est definie par une lonction non lineaire de- 
croissante, le signal d'erreur post-filtre (e2(n)) etant obtenu en appliquant une transformee de Fourier inverse aux 
composantes de Fourier ponderees (E2 p (m) ; E2 k (f)) du signal d'erreur 

12. Procede selon I'une quelconque des revendications 1 a 3, dans lequel les M composantes spectrales (z (i) (n')) du 
3$ signal d'observation sont obtenues par une decomposition en sous-bandes du signal d'observation (z(n)) t dans 

lequel la meme decomposition en sous-bandes est effectuee sur le signal de reception (x(n)) pour produire M 
composantes spectrales correspondantes (x^ (n')) du signal de reception, dans lequel le signal de prediction 
d'echo (y(n)) est produit sous la forme de ses M composantes spectrales (y (l) (n')) respectivement delivrees par 
M filtres adaptatifs auxquels sont respectivement appliquees les M composantes spectrales (x (i) (n')) du signal de 

40 reception, dans lequel le signal d'erreur (e(n)) est produit sous la forme de M composantes spectrales correspon- 

dantes (e (i) (n')) dont chacune est obtenue en soustrayant une composante spectrale (y (i) (n')) du signal de pre- 
diction d'echo d'une composante spectrale correspondante (z 0) (n'j) du signal d'observation, et dans lequel la 
grandeur (r (l) ) representative du degre de similitude entre deux composantes spectrales correspondantes (z (i) (n*), 
y (i) (n')) du signal d'observation et du signal de prediction d'echo est une valeur d'une fonction d'intercorrelation 

45 normalisee entre ces composantes spectrales. 

13. Procede selon la revendication 12, dans lequel chaque mise en forme spectrale du signal d'erreur (e(n)) comporte 
une ponderation de chacune de ses M composantes spectrales (e (i) (n')) par une quantite dependant de fagon 
non lineaire de la valeur (r (i) ) de la fonction d'intercorrelation entre les composantes spectrales correspondantes 

50 ( n *). v (i) ( n ')) du s '9 nal d'observation et du signal de prediction d'echo. 

14. Procede selon les revendications 2 et 1 3, dans lequel, pour la production du signal d'erreur filtre (e1(n)), la quantite 
de ponderation de chaque composante spectrale (e (i) (n 1 )) du signal d'erreur est une fonction non lineaire croissante 
(FNL1) de la valeur (r (i) ) de la fonction d'intercorrelation, et dans lequel chaque filtre (7) auquel est applique une 

ss composante spectrale (x (j) (n')) du signal de reception est adapte en fonction de cette composante (x (i) (n')) et de 

la composante spectrale ponderee correspondante (el^ (n')) du signal d'erreur. 

15. Procede selon les revendications 3 et 13, dans lequel, pour la production du signal d'erreur post-filtre (e2(n)), la 
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quantity de pond6ration de chaque composante spectrale (e (j) (n'» du signal d'erreur est une fonction non Iin6aire 
d6croissante (FNL2) de la valeur (r w ) de la tonction d'intercorreiation, et dans lequel le signal d'erreur post-fiftr6 
(e2(n)) est obtenu par une recombinaison, duale de la d6composition en sous-bandes, des composantes spec- 
trales pond6r6es (e2^ (n 1 )) du signal d'erreur. 

5 

16. Annuleur d'echo adaptatif, dans lequel un signal de reception (x(n)) est applique en entree d'un systeme generateur 
d'echo (1 ) et un signal d'observation (z(n)) est obtenu en sortie dudit systeme, comprenant : des moyens de filtrage 
(7) pour produire, par filtrage du signal de reception, au moins un signal (y(n)) de prediction d'un echo du signal 
de reception dans le signal d'observation ; et des moyens soustracteurs (8) pour produire un signal d'erreur (e(n)) 
w en soustrayant un signal de prediction d'echo du signal d'observation, caracteris6 en ce qu'il comprend en outre : 

des moyens (38;D') de decomposition du signal d'observation (z(n)) en M composantes spectrales (Zp(m) ; z (i) 
(n')) se rapportant k des frequences diff6rentes du spectre, M designant un nombre entier plus grand que 1 ; des 
moyens (41.m;76) de calcut de M grandeurs (msc p (m) ; <x p (m) ; r (1) ) & partir des M composantes spectrales du 
signal d'observation et de M composantes spectrales correspondantes (Y p (m) ; y (i) (n')) d'un signal de prediction 
is d'6cho (y(n)), chacune desdites M grandeurs ayant une valeur & variation progressive representative d'un degr6 

de similitude entre, respectivement, I'une des M composantes spectrales du signal d'observation et l a composante 
spectrale correspondante du signal de prediction d'6cho ; et des moyens (33,34;77,78,79,8l ,H) ae mise en forme 
spectrale du signal d'erreur (e(n)) en fonction desdites M grandeurs calcul6es. 

Annuleur d'echo adaptatif selon la revendication 16, comprenant en outre des moyens d'adaptation (9) pour regler 
dynamiquement des parametres de filtrage utilises par les moyens de filtrage (7) en fonction du signal de reception 
et du signal d'erreur, dans lequel les moyens de mise en forme spectrale comprennent un filtre (33;77,7B) auquel 
est soumis le signal d'erreur (e(n)) et dont le signal de sortie (el(n)) est fourni aux moyens d'adaptation (9). 

Annuleur d'echo adaptatif selon la revendication 16 ou 17, dans lequel les moyens de mise en forme spectrale 
comprennent un post-fiftre (34;79,81 ,R) auquel est soumis le signal d'erreur (e(n)) et dont la sortie deiivre le signal 
de sortie (e2(n)) de I'annuleur d'echo. 

Annuleur d'echo adaptatif selon I'une quelconque des revendications 16 & 18, caracterise en ce que les moyens 
(7) de filtrage du signal de reception (x(n)) sont agences pour effectuer un filtrage dans le domaine fr6quentiel. 

Annuleur d'6cho adaptatif selon I'une quelconque des revendication 16 a 18, caract6ris6 en ce que les moyens 
(7) de filtrage du signal de reception (x(n)) sont agences pour effectuer un filtrage en M sous-bandes auxquelles 
sont respectivement associees les M composantes spectrales des signaux. 
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